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SUMMARY: 

As  part  of  the  sediment  acoustics  experiment-2004  (SAX04),  bistatic  synthetic  aperture 
sonar  (SAS)  measurements  were  made.  The  Applied  Physics  Laboratory  (APL-UW) 
rail/tower  system  was  used  to  receive  signals  emitted  by  either  the  Naval  Surface  Warfare 
Center-Panama  City  (NSWC-PC)  tower  with  a  parametric  source  or  a  separate  APL-UW 
tower  with  an  omnidirectional  source.  For  the  parametric  source,  the  relevant  frequency 
bands  are  1-5  and  60-65  kHz.  The  omnidirectional  source  emitted  a  short  20  kHz  sine 
wave  burst.  The  resulting  data  has  been  preprocessed  into  matrix  form  and  is  available 
for  further  analysis  by  other  researchers. 

In  addition,  we  have  presented  the  results  of  our  analysis  to  date  of  this  same  data  in  two 
publications: 

S.  G.  Kargl,  K.  L.  Williams,  E.  I.  Thorsos,  J.  L.  Lopes,  “Bistatic  synthetic  aperture  sonar 
measurements  and  preliminary  analysis.”  in  Proceedings  of  Boundary  Influences  in  High 
Frequency,  Shallow  Water  Acoustics,  University  of  Bath,  UK  5th-9th  September  2005. 
(Copy  attached  as  part  of  this  final  report) 

K.  L.  Williams,  E.  I.  Thorsos,  S.  G.  Kargl,  D.  Tang,  J.  L.  Lopes  and  C.  L.  Nesbitt, 
“Synthetic  Aperture  Sonar  Measurements  of  Bistatic  and  Monostatic  Scattering  from 
Proud  and  Buried  Targets^  submitted  to  JUA(USN). 
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As  part  of  the  sediment  acoustics  experiment-2004  (SAX04),  bistatic  synthetic  aperture 
sonar  (SAS)  measurements  were  made.  The  Applied  Physics  Laboratory  (APL-UW) 
rail/tower  system  was  used  to  receive  signals  emitted  by  either  the  Naval  Surface  Warfare 
Center-Panama  City  (NSWC-PC)  tower  with  a  parametric  source  or  a  separate  APL-UW 
tower  with  an  omnidirectional  source.  For  the  parametric  source,  the  relevant  frequency 
bands  are  1-5  and  60-65  kHz.  The  omnidirectional  source  emitted  a  short  20  kHz  sine  w  ave 
burst  A  brief  description  of  the  measurement  geometry  will  be  given,  and  the  time-domain 
imaging  algorithm  will  be  outlined.  An  important  aspect  of  these  measurements  is  the  fixed 
location  of  the  source  towers  relative  to  the  rail/tower  system  removes  possible  problems 
with  motion  compensation  of  source  and  receiver  arrays  as  could  occur  when  the  arrays  are 
mounted  on  autonomous  underwater  vehicles.  Various  objects  were  placed  in  the  field  of 
view  of  the  SAS  system.  Imaging  results  for  one  of  these  objects  will  be  shown. 


1  Introduction 

A  bistatic  synthetic  aperture  sonar  (SAS)  scenario,  which  utilizes  multiple  autonomous 
underwater  vehicles  (AUV),  has  the  desirable  features  of  possible  rapid  high-area  coverage, 
multiple  aspect  angle  looks  at  an  object  (e.g.,  forward  scattering),  and  covert  operation. 
Unfortunately,  these  features  are  achieved  at  the  expense  of  relative  vehicle  motion. 
Traditional  towed  SAS  systems  have  implemented  motion  compensation  schemes  where  the 
source  and  receiver  are  mounted  on  the  same  vehicle.  In  the  multi- AUV,  bistatic  SAS 
scenario,  not  only  is  motion  compensation  for  each  AUV  required,  but  each  AUV  must 
register  its  position  with  all  other  AUV  positions.  This  is  required  because  coherent 
processing  of  the  recorded  signals  needs  proper  phase  information  to  form  images. 

In  this  paper,  the  problems  with  the  use  of  AUVs  are  removed,  so  that  we  can 
concentrate  on  the  multiple  aspect  angle  looks  at  an  object.  The  APL-UW  rail  system  [1] 
with  its  SAS  array  is  used  in  a  receive  mode  while  the  source  is  located  on  a  distant,  fixed, 
tower.  The  APL-UW  mobile  tower  [2]  and  NSWC-PC  parametric  sonar  tower  [3] 
transmitters  were  used  in  a  set  of  experiments  during  SAX04.  This  paper  describes  the 


measurements  and  demonstrates  that  the  resulting  pulse  compressed  baseband  signals  and 
the  bistatic  SAS  images  made  using  these  signals  yield  additional  information  for  the 
detection  and  classification  of  targets. 

2  Experimental  Configuration 

Figure  1  depicts  the  bistatic  configuration  of  the  SAX04  experiments.  Several  targets  were 
placed  on  the  seafloor  10  meters  from  the  source  tower.  For  the  measurements  with  the 
NSWC-PC  parametric  source,  the  acoustic  axis  of  the  transmitter  was  aligned  through  (or 
nearly  through)  the  proud  target  The  parametric  source  was  driven  with  two  linear 
frequency  modulated  (LFM)  chirps  where  the  primary  carrier  frequencies  were  nominally 
61  and  64  kHz  with  a  2  kHz  bandwidth  and  4  ms  duration.  One  chirp  had  a  positive  chirp 
rate  while  the  chirp  rate  of  the  second  was  negative.  The  generated  secondary  field  has  an 
approximate  frequency  range  of  1-5  kFlz.  A  second  set  of  LFM  chirps  generated  a 
secondary  field  with  a  5-16  kHz  range  where  the  carrier  frequencies  were  nominally  57  and 
67  kHz  with  a  7.5  kHz  bandwidth. 


Figure  1.  The  rail  is  located  along  the  y  axis  and  signals 
were  recorded  from  - 13.5  to  13.5  m  at  a  Ay  =  0.025  m 
interval.  This  interval  was  chosen  because  the  horizontal 
aperture  of  the  receiver  is  dh  =  0.1  m,  and  Gough  and 
Hawkins  [4]  have  suggested  that  the  interval  should  be  Ay  * 
dh/4.  The  dashed  box  denotes  the  edge  of  the  images  to  be 
shown,  and  the  largest  target  placed  in  the  center  of  the 
image  area  was  approximately  2  m  in  length. 

yW 

For  the  measurements  with  the  parametric  source,  the  NSWC-PC  tower’s  location  was  (x,y) 
=  (10,13.5)  m  with  the  source  approximately  4  m  above  the  seafloor.  The  APL-UW  mobile 
tower  had  a  nominal  location  of  (x,y)  =  (28,15)  m  and  its  omnidirectional  source  was  4.33  m 
above  the  seafloor. 

For  the  measurements  that  involved  the  NSWC-PC  parametric  source,  only  images  of 
the  target  from  the  primary  frequency  band  are  shown.  Images,  based  on  the  secondary 
field,  have  been  formed,  but  the  low  signal-to-noise  of  the  data  tends  to  lower  the  quality  of 
images.  The  strength  of  a  parametrically  generated  field  depends  on  the  amplitude  of  the 
primary  Fields.  During  SAX 04,  it  was  found  that  the  direct  arrival  of  the  primary  field  was 
sufficient  to  saturate  the  pre-amplifiers  of  the  APL-UW  receiving  array.  The  APL-UW 
electronics  were  originally  designed  for  a  series  of  experiments  that  did  not  include  the 
bistatic  SAS  measurement,  so  bandpass  filtering  prior  to  the  pre-amplification  stage  was  not 
needed.  The  addition  of  bandpass  filtering  during  the  bistatic  SAS  measurements  was  not 
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possible.  Thus,  the  amplitude  of  the  primary  fields  had  to  be  limited  to  avoid  saturation  of 
the  pre-amplifiers.  As  a  consequence,  the  strength  of  the  secondary  fields  was  also  reduced. 

3  Time  Domain  Imaging  Algorithm 

A  time-domain  image  is  constructed  via  a  delay  and  sum  algorithm.  For  a  pixel  located  at 
rij  *  (x| ,  yj ,0) ,  the  complex  image  amplitude  is 

PiJ-P(«'ij)-2A(rs,rrk,rij)p(yrlc,tk)exp[-ik0(Rs+Rrk)]  0)  where 

J  k 

t  ^  =(Rs+Rrjc)/c  is  the  time  delay  for  propagation  from  the  source  to  the  pixel  and  then 
from  the  pixel  to  the  receiver.  The  separation  distances  are  Rs  =  |  rs  - rjj  |and 

Rrk  “  lrrk  ~  rij  I  where  rs  and  rrlc  are  the  source  and  receiver  locations  and  ry  is  the 
pixel  location  on  the  mean  sediment  surface  z  =  0,  and  the  speed  of  sound  in  the  water  is  c  = 
1538  m/s.  The  amplitude  factor,  ACr^r^ry),  contains  the  beam  patterns  of  the  source 

and  receiver  as  well  as  a  time  varying  gain  to  enhance  the  contrast  across  an  image.  In  the 
formation  of  the  images  in  Sec.  4,  contributions  that  occur  outside  of  the  -6dB  down 
points  of  the  main  lobe  of  the  beam  patterns  are  ignored  (i.e.,  suppression  of  possible 
interference  from  side  lobes).  The  complex,  pulse-compressed,  baseband  pressure  at  the 

receiver  location  is  p(yr|c,tjc).  Phase  factor  in  Eq.  (1),  exp[-ik()(Rs  +Rrjc)],  is  phase 

compensation  for  a  spherically  diverging  wave  and  the  assumed  time  convention  of 
exp(-io)t)  where  co  and  k0  are  the  angular  carrier  frequency  and  wavenumber  in  water, 

respectively.  If  the  y-axis  is  aligned  with  the  along-track  direction,  then  the  separation 
distances  are 

Rs  -  tx? + (y j  -  ys)2  + zs  ],/2>  Rrk  =[*?  +(yj-yrk)2+zr2]l/2.  (2>  where 

the  image  plane  is  tacitly  assumed  to  be  at  z  =  0;  however,  this  restriction  is  not  necessary 
and  other  image  planes  can  be  defined  by  the  obvious  change  to  Eq.  (2).  It  is  also  noted  that 
Eq.  (2)  implies  that  only  the  receiver  moves  along  the  SAS  track  via  the  k  subscript.  To 
form  an  image,  the  discrete  representation  of  the  spatial  coordinates  are 

X,  “(x0 -xL/2)  +  (i-l)Ax,  (i  =  l,---,Nx)  (3) 

yj-(yo-yL/2)+(j-l)Ay,  (j-l,-,Ny)  (4) 

y* -yrO+(k-])Ayr,  (k  =  l,--,Nr)  (5)  where 

the  center  of  the  image  is  at  (xo,y0),  the  length  of  the  sides  of  the  image  are  Xl  and  yL,  and 
the  initial  receiver  location  is  y^. 

Inspection  of  Eqs.  ( \)-(5 )  and  the  discrete  nature  of  the  sampled  time  signals 
p(y*>  tn)>  (n  =  I,...,  N,)  suggests  that  t*  seldom  coincides  with  a  discrete  time  point.  Thus, 
the  summation  in  Eq.  (1)  requires  interpolation  of  the  signals.  The  interpolation  has  been 
performed  with  6,  12,  24,  and  48  point  cubic  splines,  which  were  constructed  on  a  discrete 


time  interval  containing  tk.  Although  cubic  spline  interpolation  permitted  the  construction 
of  an  image,  a  4  point  Lagrange  interpolation  algorithm  was  found  to  be  sufficient  for  image 
formation  and  the  Lagrange  interpolation  is  more  numerically  efficient  than  cubic  spline 
interpolation.  For  the  results  presented  here,  Lagrange  interpolation  was  used. 
Additionally,  the  recorded  signals  were  sampled  at  1  MHz,  and  the  use  of  the  pulse- 
compressed  baseband  representation  yields  a  finely  sampled  signal,  which  significantly 
reduces  the  potential  for  numerical  round-off  errors  in  computing  the  interpolated  complex 
pressure. 


4  Discussion 

Figure  2  shows  the  magnitude  of  the  pulse-compressed,  baseband  scattered  pressure  and  the 
image  formed  via  the  time  domain  algorithm  described  above.  The  target  is  an  open-ended 
cement  pipe,  which  is  nominally  1.1  m  in  length  and  0.42  m  in  diameter.  The  data  were 
taken  on  21  October  2005.  In  Fig.  2 A,  the  rail  system  is  located  at  the  top  of  the  image  and 
the  parametric  source  is  to  the  right  of  the  image.  The  earliest  arrival  that  forms  a  nearly 
parabolic  arc  is  the  direct  arrival  from  the  source  to  the  receiver.  The  next  parabolic  arc  is 
the  first  bottom  bounce  of  the  source.  The  bright,  broad  arc  that  appears  after  13  ms  is  the 
scattering  from  the  cement  pipe.  It  is  noted  that  the  y  =  13.5  and  —13.5  m  locations  are  at 
approximately  45°  and  150°  with  respect  to  backscattering  to  the  source.  As  the  receiver 
moves  into  the  near  forward  scattering  region,  the  direct  arrival,  bottom  bounce,  and 
scattered  signal  nearly  simultaneously 
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Figure  2.  The  orientation  of  the  cement  pipe  gives  an  end-on  ensonification.  These  images  were 
formed  from  the  LFM  chirps  of  the  primary  frequencies  and  each  image  is  normalized  to  its 
“hottest”  pixel.  A.  Pulse-compressed  baseband  signals.  B.  Bistatic  SAS  image.  The  black 
rectangle  represents  the  nominal  dimensions  and  orientation  of  the  cement  pipe. 


Figure  3.  The  cement  pipe  is  oriented  in  a  broadside  configuration.  A.  Pulse-compressed 
baseband  signals.  B.  Bistatic  SAS  image.  The  black  rectangle  represents  the  nominal 
dimensions  and  orientation  of  the  pipe. 

reach  the  receiver.  Figure  2B  shows  the  image  formed  from  the  data  in  Fig.  2A.  Note,  the 
rail  system  is  located  at  the  bottom  of  Fig.  2B  One  can  identify  the  front  and  back  ends  of 
the  pipe  where  the  bright  feature  on  the  right  is  nearest  the  parametric  source.  The  length 
and  width  of  the  pipe  can  be  estimated  from  this  image  to  be  on  the  order  of  1  m  long  and 
0.5  m  wide.  The  spread  in  the  image  may  be  due  to  multiple  scattering  of  acoustic  energy 
within  the  open  ended  pipe. 

Figure  3  shows  the  magnitude  of  the  complex,  pulse-compressed,  baseband  pressure 
and  the  image  formed  from  the  data.  The  data  were  recorded  on  22  October  2005.  The  ends 
of  the  pipe  are  again  evident  in  Fig.  3B  at  ranges  of  10  and  1 1  m,  respectively.  However, 
the  brightest  feature  appears  at  the  center  of  the  pipe.  The  width  of  the  pipe  appears  smaller 
than  0.25  m,  which  is  partly  due  to  the  fact  that  neither  end  is  directly  illuminated  with 
sound.  As  in  Fig.  2A,  Fig.  3A  shows  the  direct  arrival,  bottom  bounce,  and  the  scattered 
acoustic  signals. 

The  NSWC-PC  parametric  source  was  designed  to  have  a  narrow  beamwidth  for  its 
secondaiy  field,  and  hence,  the  primary  field  also  has  a  narrow  beamwidth  (i.e.,  on  the  order 
of  5°)  [5].  In  a  bistatic  SAS  scenario,  it  may  be  tempting  to  light  up  the  underwater 
environment  with  an  omnidirectional  source  to  obtain  high-area  coverage.  This  was  carried 
out  with  the  APL-UW  mobile  tower  on  16  October  2005.  A  three-cycle  20  kHz  sine  wave 
pulse  was  broadcast  from  an  omnidirectional  transducer  on  the  mobile  tower  with  the 
approximate  position  of  (28,15,4.33)  m. 

Figure  4  depicts  the  magnitude  of  the  complex,  pulse-compressed,  baseband  pressure 
and  the  image  formed  from  the  data.  The  receiver  was  moved  along  the  rail  at  an  of  interval 
Ayr  =  0.05  m.  The  increased  interval  was  due  to  limitations  in  the  signal  generation  for  the 
mobile  tower  source.  The  baseband  pressure  contains  four  broad  parabolic  features  that 
span  the  cross  range  dimension.  These  are  the  direct  arrival,  the  bottom  bounce,  the  surface 
bounce,  and  the  bottom-surface  bounce.  With  a  nominal  water  depth  of  18  m  and  receiver 
location  of  (0,13.5,3.91)  m,  these  arrivals  should  appear  near  18.2,  19.0,  25.7,  and  30.7  ms, 
respectively.  The  slight  discrepancy  of  the  computed  arrival  times  and  those  observed  in 
Fig.  4 A  can  be  attributed  to  uncertainties  in  mobile  tower  location  and  water  depth. 
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Figure  4.  A.  Pulse-compressed  baseband  signals. 


B.  Bistatic  SAS  image  with  three  targets. 
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The  important  features  in  Fig.  4A  are  the  thin,  crescent-shaped  arrivals  near  25  ms  and 
cross  range  of  5  to  13.5  m  and  near  28  ms  and  cross  range  of  2  to  13  m.  These  indicate 
objects  are  present,  and  Fig.  4B  reveals  at  least  two  targets  near  (4,31)  m  and  (8,31)  m.  A 
third  target  may  be  located  near  (9,30)  m.  These  targets  coincide  with  several  proud  targets 
placed  on  the  seafloor  for  other  SAS  measurements 


5  Conclusions 

The  combination  of  a  bistatic  SAS  image  and  an  image  of  its  pulse-compressed  baseband 
data  may  provide  additional  information  for  classification  and  identification  of  a  target  that 
may  not  be  present  in  data  from  conventional  monostatic  SAS  systems.  Targets  with 
complex  shapes  may  have  low  backscattering  target  strengths,  but  relatively  high  bistatic, 
and  in  particular,  forward  scattering  cross  sections.  Inspection  of  Figs.  2A  and  3A  shows 
strong  scattering  from  the  cement  pipe  in  the  angular  range  of  approximately  90°  to  145° 
with  respect  to  backscattering.  The  high  signal-to-noise  ratio  provides  the  possibility  of 
processing  schemes  to  extract  elastic  responses  of  the  target.  Hence,  alternative  processing 
schemes  may  be  able  to  extract  this  information. 

Bistatic  SAS  images  were  constructed  from  data  obtained  during  SAX04.  The  images 
in  Figs.  2B  and  3B  clearly  indicate  the  presences  of  a  target  and  an  estimate  of  its  size  can 
be  achieved.  Even  though  these  images  are  substantially  different  due  to  the  orientation  of 
the  cement  pipe  with  respect  to  the  source  location  and  rail  system,  unambiguous 
identification  of  the  cement  pipe  from  these  bistatic  SAS  images  is  unlikely.  Bistatic  SAS 
images  alone  may  offer  advantages  over  conventional  monostatic  SAS  images  for  detection 
due  to  higher  bistatic  target  strengths.  Also,  after  a  target  has  been  detected,  the  same 
bistatic  data  can  give  important  information  on  classification.  Figures  2A  and  3 A  show  the 
magnitudes  of  the  complex,  pulse-compressed,  baseband  pressure.  The  parabolic  arcs 
associated  with  scattering  from  the  target  in  the  two  orientations  are  clearly  different,  and 


these  may  contain  target  information  (e.g,,  elastic  resonances)  that  is  important  to 
classification  and  identification  of  the  target.  SAS  processing  removes  this  information. 

The  complex,  pulse-compressed,  baseband  scattered  pressure  shown  in  Fig.  4A  for  an 
omnidirectional  source  clearly  demonstrates  the  need  to  resolve  weak  target  scattering  from 
multipath  arrivals.  Otherwise,  the  signal-to-noise  ratio  may  be  too  low  to  adequately  resolve 
a  target.  For  example,  a  third  target  may  be  located  at  (9,30)  m  in  Fig.  4B.  One  strategy  to 
achieve  high-area  coverage  would  be  to  have  broad  horizontal  and  narrow  vertical  beams 
directed  at  a  small  depression  angle  to  the  surface. 
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